The general Penrose inequality: lessons from numerical evidence. 
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Formulation of the Penrose inequality becomes ambiguous when the past and future apparent 
horizons do cross. We test numerically several natural possibilities of stating the inequality in 
punctured and boosted single- and double- black holes, in a Dain-Friedrich class of initial data and 
in conformally flat spheroidal data. The Penrose inequality holds true in vacuum configurations for 
the outermost element amongst the set of disjoint future and past apparent horizons (as expected) 
and (unexpectedly) for each of the outermost past and future apparent horizons, whenever these 
two bifurcate from an outermost minimal surface, regardless of whether they intersect or remain 
disjoint. In systems with matter the conjecture breaks down only if matter does not obey the 
dominant energy condition. 

Dedicated to Andrzej Staruszkiewicz on the occasion of his 65"" birthday. 



I. INTRODUCTION. 

The Penrose- Hawking (P], 0) singularity theorems point at the incompleteness of the classical general relativity. 
The cosmic censorship hypothesis can be regarded as an attempt to contain the damage, by demanding that the 
genuine singularities are hidden within the black holes. While it is still not clear whether the cosmic censorship 
hypothesis holds true, there is no doubt that it has shaped the research field and led to many important results 
concerning evolving systems (for a recent review see Penrose invented in 1973 an inequality that can constitute 
a necessary condition for the validity of the cosmic censorship (j^, [^). The Penrose inequality has been recently 
proved, following a scenario suggested in 1973 by Geroch in the important riemannian case (0, 0), that can be 
roughly described as a momentarily static data set of Einstein evolution equations. It was known for a long time to 
hold in spherically symmetric systems [lOj ; a later analysis allowed one to elucidate the problem of the needed energy 
conditions and show that it is independent of the foliation conditions (^3, There exist analytical scenarios 

for the proof of the general Penrose inequality (0| and but their validity is not proven, due to the technical 
complexity. 

While numerics cannot per se produce a proof, it may disprove the hypothesis or, more likely, be of help in finding 
its correct formulations. The latter is the main goal of this paper - Sec. 3 presents possible wordings of the Penrose 
inequality that are checked in later sections in examples representing several classes of initial data. The obtained 
results can be described as bringing some surprises. In the spherically symmetric spacetimes the future and past 
apparent horizons are disjoint and it is well known that the Penrose inequality is valid only for the outermost of all 
apparent horizons (later on abbreviated as AH) assuming an energy condition. In nonspherical spacetimes the 
two AH's can intersect in most space-time foliations, including the maximal ones. Notable exceptions are the polar 
gauge foliations in which the two horizons cannot be separated (although they can bifurcate), but their existence status 
is unclear. The apparent ambiguity would be resolved by accepting the proposal of Horowitz [l5l | that one should 
take a surface of a minimal area enclosing all AH's. It is unexpected in this context that whenever the outermost 
past and future apparent horizons bifurcate (in a sense specified in Sec. 5) from an outermost minimal surface, then 
the inequality holds true for each of them. This remark applies to either crossing or disjoint AH's, in all nonspherical 
vacuum configurations tested by us. One can convert this phenomenological observation into a local analytic proof, 
as sketched in Sec. 5. 

The order of remainder of this paper is following. Next section brings the initial constraint equations and a brief 
description of the conformal method of constructing initial data. Several versions of the Penrose inequality are given 
in Sec. 3. Sec. 4 briefly describes the relevant numerical methods. Obtained results are reported in Sees 5-7. They 
support the various versions of the inequality for the vacuum initial data in Sees 5 (punctured-Bowen-York data) and 
6 (punctured Dain-Friedrich data). Section 7 deals with non- vacuum spheroidal initial data; in this case the Penrose 
inequality can be broken, if the dominant energy condition is not valid. Last section presents main conclusions. 



II. EINSTEIN CONSTRAINT EQUATIONS. 



Let S be an asymptotically flat Cauchy hypersurface endowed with an internal metric gij, the scalar curvature R 
and the extrinsic curvature Kij. The initial constraint equations read |0| 

R = IQup + KijK'^ - {Kff 

W4K}-g}Kl) ^Sttj,, (1) 

where p and j are the mass density and current density of initial material fields. In the case of maximal slicing 
condition, = 0, the initial data can be found by the conformal method 17]. In what follows we analyze conformally 
flat classes of solutions, corresponding to vacuum and spheroidal systems with matter. The metric reads gij = (jy^gij 
where gij is the Euclidean metric. 

The global energy-momentum can be found from standard formulae 



E ^ — I SS'W, 



-1 

1 



Pi = izl d^s^K^v (2) 



The asymptotic mass is given by m = ^E'^ — PiPK 

An apparent horizon will understood later on as a two-dimensional surface S lying in S with a normal t satisfying 
one of the two equations 

e+ = V^f^K.jtH^ = 0, (3) 

where the sign and _ corresponds to the past and future apparent horizons, respectively and 0's are known as 
optical scalars. 



III. FORMULATION OF THE PENROSE INEQUALITY. 

The Penrose inequality is expected to hold only for the outermost AH. As explained in in the case of spherical 
symmetry: Consider the outermost future trapped surface, the (future) apparent horizon, call it S. Let us assume that 
S is outside the outermost past trapped surface. In other words, we assume 0+{S) = and that both 6^ and 0_ are 
positive outside S. A simple analytic argument shows the validity of the following inequality 

provided that the dominant energy condition is satisfied by matter located outside S. Sr in this formula is the area 
of S. As stressed in the quoted paper. Of course, an identical argument works if the outermost trapped surface is a 
past apparent horizon. 

In the spherically symmetric geometries optical scalars have the same level sets, since both of them are constant on 
centered spheres. Thus the spheres surrounding the outermost AH have positive 9- and 0+ and i) possess larger area 
than the AH. Obviously ii) the future and past horizons do not cross. Therefore there is no ambiguity in defining the 
Penrose inequality and - since the cosmic censorship hypothesis asserts that AH's are enclosed by the event horizon 
that asymptotically evolves to a Schwarzschild or Reissner-Nordstroem black hole horizon - it can be regarded as the 
necessary condition for the cosmic censorship ^ . 

None of the features i) and ii) becomes obvious in the general nonspherical case, even if our liberal definition of 
AH's is replaced by the more stringent one due to Penrose. Let us recall that in Q future trapped surface are assumed 
to have - in our terminology - a positive scalar 0_ and a negative scalar 9^. That is, each of the two beams of null 
geodesies emanating orthogonally outward and inward from a trapped surface, is convergent. Consequently, a future 
apparent horizon (understood as the outermost boundary in the set of all future trapped surfaces) has vanishing 9+ 
but non-negative 9-. The analogous situation (but with optical scalars reversing their roles) takes place for past 
apparent horizons. Thus the picture of AH's that emerges here resembles that of spherically symmetric geometries. 
There exist level sets of, say 9+ < 0, such that 9- > Q (and conversely). The two optical scalars do not possess 
common level sets, but the sign of one of them is controlled on the level set of the other. This can happen only if one 



matches in a suitable way the choice of both a Cauchy hypersurface and of the two-dimensional foliation within this 
slice. (This is inherent also to the scheme of the proof of the Penrose inequality that is proposed in One can 

find a two-surface S to be such an AH in one particular foliation, but that may not be true in other space-like slices. 
Even with this stringent definition, the future and past AH's can intersect and there may exist surfaces of a smaller 
area surrounding them as pointed out by Horowitz [l5j . 

In the rest of this paper by AH's are understood two-surfaces satisfying one of the conditions of Q, which might 
be weaker than the notion employed in the singularity theorems (but see a discussion following the point ii) below). 
That means that the failure of a particular version (or all of them) formulated below of the Penrose inequality does 
not necessarily negate the the cosmic censorship hypothesis (CCH). And conversely, their validity lends even more 
credence in CCH. It is not without significance that such AH's are easier to find numerically than the standard objects 
defined by Penrose. 

The three versions of the Penrose inequality read as follows, (assuming the dominant energy condition for 
nonvacuum initial data) 

i) The minimalistic one (PIM henceforth); it was borrowed from a proposition first put forward by Horowitz 
The surface Am of a smallest area Sm surrounding regions with horizons satisfies the inequality 



It appears in many of the numerical cases reported later that Am coincided with that constructed from apparent 
horizons (see ii)), but in a number of configurations it consisted also of segments of minimal surfaces. The existence 
of configurations having portions of minimal surfaces extending outside AH's, means that it is not excluded that the 
actual area of an event horizon - if there is one - is smaller than that of the AH. In such a case PIM constitutes the 
necessary condition for the validity of CCH. 

ii) The standard one (PIS later on). The closed 2-surface Ah is either the outermost apparent horizon (if AH's do 
not intersect) or a union of segments of outermost future and/or past apparent horizons. Then its area Sh satisfies 
the inequality 



Ah does not manifestly satisfy the assumptions of the singularity theorems, but its importance lies in the fact that 
it may do so in another foliation (say, polar gauge one). The heuristic argument is as follows. The product of two 
optical scalars is i) invariant and ii) vanishes on Ah- If there exists a local boost to a polar gauge foliation (that is, 
a foliation with 9- —9+] apparent horizons correspond here to minimal surfaces) of a space-time, then on a polar 
gauge slice the two-surface Ah would become an apparent horizon (that is, the outermost minimal surface) in the 
sense of Penrose, and then the CCH demands the existence of an event horizon. The area of the intersection of the 
event horizon with the actual polar gauge slice would have to be bigger than of Ah- Accepting that, PIS seems to be 
just right one, as a necessary condition, from the point of view of CCH. Unfortunately, there is a gap in the argument. 
Namely, there is no possibility to rule out the existence of minimal surfaces that extend outward of the outermost 
apparent horizon. There are reasons to expect (basing on the analogy to spherical symmetry) that the polar gauge 
slice does not penetrate regions with minimal surfaces and the surface Ah would not be seen on the slice. 

iii) In the cases with intersecting apparent horizons Aa's we compared also their area related quantities \J~^^ with 
the asymptotic mass. Invariably it was found that 



for each of the horizons, and with a significant safety margin. There is no obvious reason why the area Sa should be 
bigger than Shi but this is what was found to be true in all analyzed examples. That suggests that iii) is stronger 
than PIS and it comes as a surprise that numerics supports the inequality (|7|. 

Data concerning the two stronger (ii) and iii)) of the above conjectures are given in sections 5 and 6; their validity 
implies PIM. Only in the first part of Sec. 5 we report data concerning the version PIM, to show that the results are 
close for all three statements. The horizons do not cross in the case of spheroidal initial data and there is no minimal 
surface outside the outermost AH; therefore conjectures PIS and PIM do coincide in the examples considered in Sec. 
7. 



IV. DESCRIPTION OF NUMERICAL METHODS. 



In the conformal method and for vacuum conformaUy flat initial data, one first solves analytically the equation 
WiKj = (the covariant derivatives are in the Euclidean metric) and then the Lichnerowicz-York equation 

A0^-:^0-^ (8) 

with the flat laplacian A. This is an example of a weakly nonlinear elliptic equation; its leading derivatives are linear 
(hence the equation is quasilinear) and the nonlinearity is rather weak (c.f. negative powers of the conformal factor 
(/)). Due to the cylindrical symmetry, we search for as a function of the angle 6 and the coordinate radius r. It 
is necessary to map the problem into one with finite domain; thus the radius r is replaced by another independent 
variable v = r/(l + r). Adopting x = cos 9, one has to solve (jSj) in the rectangular — l<a;<l,0<t'<l. It is solved 
iteratively by the standard Newton method on the lattice up to 200 x 5000 points. Due to the weak nonlinearity 
of the Lichnerowicz-York equation, it was enough to apply at most 4-5 iterations. We used four different solvers, in 
particular the MUMPS 1^ and HYPRE ^ ones. 

The apparent horizon equation (O becomes in our context a nonlinear ordinary equation, for the function r(9). This 
is a classical two-point problem (see a discussion in |23| in a similar context) with dr / d6\g=Q = dr/d9\e=T^ = 0. It is 
solved by the standard shooting method. We resorted to two numerical packages, ODEPACK |2ll| and SUBPLEX l2^ . 
The needed extrapolation of the formerly found solution cf) has been done with the help of the bilinear interpolation 
|23| . The bilinear method appeared entirely satisfactory, due to the high density of our numerical lattice. 



V. BOOSTED PUNCTURED DATA. 

Assume P to be a constant vector and n - a unit normal to a metric sphere in the Euclidean geometry. Assume 
standard spherical coordinates r, 9 and (p. One can easily check that the extrinsic curvature 

satisfies the momentum part of Eq. with the vanishing current (that is the boosting part of the Bowen York initial 
extrinsic curvature (24] ) . The hamiltonian constraint (the first equation in Q) reads now, assuming vacuum case 
(p = 0) and aligning the z-axis along P, 

A^ = -^(l + 2cos20)r^ (10) 

where A is the flat laplacian. For these boosted data one obtains the global momentum Pi = Pi, the global energy E 
is given by (01 and the asymptotic mass reads ^/ E"^ — P"^ . 

There are two established ways of solving the resulting (Lichnerowicz-York) equation, i) In the first approach, that 
takes care about the global topology of the manifold (the so-called conformal imaging method 24]; that actually 
requires the use of a larger set of extrinsic curvature data), Bowen and York solve Eq. H10() outside r > a, assuming 
that the sphere r = a is a minimal surface and that at infinity the conformal factor goes to 1. ii) In the second 
approach, the puncture method, one splits into two parts. (}ii = 1 -I- mi/(2r) -f 0i and finds a solution 0i in the whole 
Euclidean space, demanding that at infinity <p\ w di/(2r) l25|; here mi and d\ are some constants, that are related 
to the global energy of the manifold. 

In this paper we use the second approach in order to treat vacuum configurations with one or two black holes. The 
data for three exemplary configurations (chosen from a much bigger sample) with a single black holes are presented 
in the first table. The first column is the parameter m,\ appearing in the preceding splitting, the second column is the 
linear momentum P, and the third, fourth and fifth columns, respectively are the global mass m and the "horizon" 

masses m^/ = ^ ran = \J^^ and niA = \J j^- In the case of single boosted black holes the asymmetry causes 
the horizons to intersect the minimal surface, and this is why we consider the case PIM with Am- The surfaces in 
the fourth and fifth columns do not coincide with apparent horizons, but consist of two {Ah) parts (of an apparent 
horizon to the future and to the past) or three segments {Am) (of an apparent horizon to the future, a minimal 
surface and an apparent horizon to the past) - as explained in Sect. III. It happens that the areas Sa of the apparent 
horizons to the past and to the future are equal; the sixth column brings corresponding data which (unexpectedly) 
obey the Penrose inequality. Each row describes a different configuration. 
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8.061855 


8.059402 


8.059426 


8.05948097 
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4.122407 


4.110666 


4.110751 


4.11093771 


4 


5 


4.707092 


4.499335 


4.500002 


4.50143757 



In this case areas of the past and future AH's are equal. It is clearly seen that all versions, the weaker (PIM) and 
the stronger (PIS) as well the last one (iii) of the Penrose inequality are satisfied. It is noticeable that the areas of 
Am^ Ah and A a are very close. The fact of interest is the numerical evidence for the existence of parts of minimal 
surfaces that extend outward of outmost apparent horizons. 

The corresponding results for two black hole configurations are comprised in the next table. Now the puncture 



method requires that (j) ~ 1 
curvature reads 



"^j/(2|7^— ri\) + (j)i (the two black holes are located at ri, i=l,2). The extrinsic 



2r202 



E 



n ) 



(11) 



where P^^'^ aligned along the 2;-axis can be interpreted as the linear momentum of the s — th black hole and n'*^ — 
(r — Rs) /\r — Rs\. As before, one finds a solution 4>i of the Lichnerowicz-York equation in the whole Euclidean space, 
demanding that at infinity (pi ~ di/(2r) 25j; as before mi and di are some constants, that are related to the global 
energy of the manifold. The first and second columns describe parameters ("mass" mi and "momentum" P*^^^ of the 
first black hole, the third and fourth columns give the same information about the second black hole. The fifth and 

sixth columns, respectively are the global mass and the areal mass mn — \J j^- Sh is the area of Ah, the 2-surface 

constructed according to the recipe ii) of Sec. II. The 2-surface Ah surrounds both black holes (located at r = 1 
and 6* = or 6* = tt). In two cases (fifth and ninth) the past and future AH's do cross; the seventh column presents 

relevant values of niA — \J j^, where Sa is the larger of the two areas in question. They satisfy all formulations of 
the inequality. 
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0.25 


10.000103 


9.999632 
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10.010290 


10.008087 
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9.006743 


9.005594 
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-2.5 
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0.5 


10.052687 


10.051593 


10.051593 



In the remaining seven cases the past and future horizons do not cross and there is a minimal surface in between 
them. Surprisingly - and in a sharp contrast with the corresponding case in spherically symmetric configurations - 
the Penrose inequality is valid simultaneously for past and future AH's. Another interesting observation is that, in all 
cases, when the parameter P tends to zero then AH's tend to the minimal surface. In this sense, the AH's bifurcate 
from the minimal surface. This is true, as matter of fact, in all numerical examples studied in Sees. 5 and 6. It 
happens, that there always exists at least one minimal surface; those AH's horizons that bifurcate from the outermost 
one do satisfy the inequality. On the other hand, AH's that branch from an innermost minimal surfaces (there are 
several such cases in our sample of data) do break all aforementioned versions. 

We show below an analytic argument that these observations remain true (with some reservations) for initial data 
with AH's that arise from small perturbations of data with minimal surfaces. 

Theorem. Let I be a real parameter, IKij be the Bowen-York or Dain-Friedrich extrinsic curvature (multi- 
puncturized) and 4>i be a solution of the Lichnerowicz-York equation on 



(12) 



Then (pi, K^j = Kij/cpf constitute initial data of the Eintein equations; assume that for each I there exist apparent 
horizons that in the limit I — > coincide with a nonspherical outermost minimal surface Sq. Then there exists Iq such 
that for |Z| < Iq the Penrose inequality is satisfied. 



Sketch of the proof. In the case of a nonspherical surface So given by r — ro{9), one has a strict inequality, 
e = niQ — yMfl7(16vr) > 0. One can easily show that: i) (using arguments of 17]) (pi > (po and mi = toq + ciP; 
ii) (using the Green function of the flat laplacian) pi < tpo + ■ Here and below Ci {i — 1, 2, ...) are some positive 
constants. The apparent horizon equations 0+ = or = depend on I through the extrinsic curvature terms 
and the conformal factor p. Since by assumption horizons bifurcate from So, they must be located within annulus 
{ro{0,(p) — czl,ro{0,(p) + C4Z). (This is due to the implicit function theorem). At an apparent horizon one should 
compare to/ = toq + ciP with the area of AH's, which is bounded from above by Ao + C5Z. It is clear that by choosing 

^0 small enough one can ensure that to/ > \J j^. 

This proof is insensitive on the sign of / and therefore it is valid for both past and future apparent horizons. Any 
attempt to convert this local result into global one would have to be preceded by a careful estimate of the dependence 
of the location r{9) of an AH on the bifurcation parameter I. Notice that this theorem does not apply to single- 
puncture initial data, since in this case the geometry corresponding to ^ = is spherically symmetric and e = 0. On 
the other hand, this result should hold for the two-puncture solutions. There is also a possibility of generalizing the 
above onto case with initial data given in the exterior of a two-surface 5i, instead of . 



VI. DAIN - FRIEDRICH CONFORMALLY FLAT INITIAL DATA. 



The main feature of these initial data is that the spatial part of the metric is conformally flat (as before), the 
momentum flow density vanishes (again, as before) and the extrinsic curvature is given, in spherical coordinates, by 
|the forthcoming formulae are translated from the language of the Newman-Penrose formalism, originally used in 



Kf = 0, for i 1,2 
^ dlW 



r^pe ^ 
1 

sin 9r^p^ 
1 

sin'^ 9d 



drd^W 



dr (rdrW) + - {xd,^W - W) 
r 



(13) 



Here x = cos 9 and W is an arbitrary function of r and 9. The component can be found from the maximal slicing 
condition Kl = 0. Let us point out that this solution generalizes the Bowen-York solution of momentum constraint; the 
latter corresponds to a particular choice of W 27]. The data of for instance, correspond ioW = rP {x^ — 3xW2. 
Extrinsic curvature (|13|l constitutes a partial case of the general solution found by Dain and Friedrich in 2001 26]. 

The conformal factor (p satisfies the Lichnerowicz-York equation = —^KijK^^p^. We seek a solution 0, using 
the puncture method, of the form p = 1 + TOi/(2r) + pi. 

The numerical calculations have been performed in the following cases: 

i) dxW — Pr {—X + x^). It is noticeable that here the global momentum is nonzero. Numerical results are given in 
the forthcoming Table. 



TOi 


P 


TO 


rriH 




4 


5 


5.622456 


4.296457 


4.296680 
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4.078569 


4.015156 


4.015260 



ii) d^W = Pr{l 
results. 



x'^)(3x'^ — l)/2. In this case the global momentum vanishes. The Table presents the obtained 
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6.073919 


4.366950 


4.372253 
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4.109128 


4.021109 


4.021207 



In both cases i) and ii) the two stronger versions, PIS and iii), of the Penrose inequality holds true. The 2-surface 
Ah, whose areal mass is depictured in the last but one column, is built from many sections of the intersecting past 
and future horizons; the number of the intersections seems to depend (for a given nonzero P) on the shape of W as 
a function of 9. There are more crossings in case i) (three ) than in the case ii) (only two). 



VII. SPHEROIDAL SYSTEMS WITH MATTER. 



Assume a foliation of the Euclidean space by oblate spheroids, 

,2 ^2 



^ ■ :^ = i- (14) 



The variable a changes from to cx3, and angle variables are r (changing from -1 to 1) and </> (varies as usual from 
to 27r). Assume that there exists a normal flow of matter with the only nonzero component 

Let hi denote the unit normal (in the Euclidean metric) to a spheroid. The related traceless extrinsic curvature reads 

At^-r»(.w49;)^^j^; (16) 

this pair, Kij and ji, solves the momentum constraint part of Eq. Let us remark, that one could always add 

to the extrinsic curvature the diagonal components K!^ — C/ {(f)^r^), Kg — K'^ — (where C is a constant and 

r = ^yx^ + y^+z^) without changing the momentum flow. We do not do this, because our primary intention is to 
study the influence of the energy conditions onto the validity of the Penrose inequality, and the aforementioned part 
of the extrinsic curvature is irrelevant from this point of view. 

The energy density p can be chosen in an arbitrary way, but the simplest possibility - that eases the analysis of the 
energy conditions - is to assume 

P = C.^-^.^ (17) 

87^(^72 + l)(a2+r2)5/2 

Later we shall put either C = 1 - which ensures the dominant energy condition - or C = 0, which breaks the energy 
conditions. Notice that 

K,iK'' = ^ X (j)-^^ X ^. (18) 

The Lichnerowicz-York equation takes now the form 



A0 = - 



1 



C 



1 



12 (a2 + T^f 



.5/2 



(19) 



4 (a2 + 1) (^2 , , , 

This equation has been solved with (p tending to 1 at infinity and bearing a constant value at the inner boundary, that 
is assumed to be a unit sphere r = 1 in the background Euclidean geometry. In the examples shown below, horizons 
do not intersect and in all cases the outermost apparent horizon (future or past) has been located outside the minimal 
surface. (There appear also initial data with intersecting horizons, but they are not of particular interest.) Its area 

Sh enters the forthcoming data through the formula mn ^ 
We shall present data corresponding to: 
i) C=l, 



0(r = 1) 
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rriH 
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2.801054 


2.800734 


2.5 


3.0009593 


3.0006660 


3.5 


5.0004402 


5,0003797834. 



It is clear that itih > rn and that the version PIS holds true. 

ii) C = 0. The energy density vanishes and therefore the dominant energy condition is broken. One expects that the 
Penrose inequality (PIS) may be broken now, and in fact this is what happens, albeit only in the first two examples. 
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2.80004295 


2.80011252 
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3.00003632 


3.000604868 


3.5 


5.000009651 


5.000005477 



VIII. CONCLUDING REMARKS. 



The weakest form of the Penrose inequality due to Horowitz is that m > -^/jm/IBtt, where Sm is the smahest area 
of a two-surface encompassing a region with apparent horizons satisfying the assumptions of the singularity theorems. 
This paper deals with three other formulations, of which even the weakest (PIM) is stronger than the Horowitz's 
one, because our notion of the outermost apparent horizon is weaker than that required by the singularity theorems. 
Despite this fact, all investigated statements of the Penrose inequality are confirmed by our numerical analysis for 
vacuum initial data and for those systems with matter that satisfy an energy condition. The only negative examples 
correspond to data with matter that does not satisfy an energy condition. 

There exist minimal surfaces in all investigated examples with vacuum; these are mostly singlets but in a number 
of cases also doublets. It is observed that the two (past and future) AH's, that bifurcate (with the momentum being 
the bifurcation parameter - see the end of Sec. 5) from the outermost minimal surface, do satisfy - regardless of 
whether they cross or do not cross - all versions of the Penrose inequality. An analytic proof that this is true (with 
some reservations), at least for initial data with small extrinsic curvature, is sketched at the end of Sec. 5. 

In summary, results of this work show that the Penrose inequality does not hold only when expected not to hold; 
that points strongly in favour of the validity of the Penrose's conjecture in physically interesting cases. 
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